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ABSTRACT 


The present work on substitution deals with various 
models where the interaction of demand due to shortages and/or 
introduction of new items is taken into account in two— item/ 
multi-item inventory system. 

Deterministic and stochastic models have he® developed 
for various situations. In deterministic case* two-way complete 
substitution of two items of decaying nature is considered 
and. a solution methodology is suggested. A simulation model 
to solve the problem of single it® one-way substitution by 
n ain&ainiag end-item substitute and second level item inventories 
as safety stocks along with emergency ordering: .for end-item 
substitute to meet a particular service level is developed. 

A simulation model for two-way substitution of N-items with 
periodic review policy is studied and an iterative method to 
find the decision variables of the model is described. 

Solution methodology has been suggested in all cases, 
aqne* 

Numerical examples/demonst rated to show the steps involved. 



CHAPTER I 


INTRODUCTION 

1.1 SUBSTITUTION - AN OVER VIEWS 

Organisations try to boost sales by improving quality 
of the product at a lower price and by maintaining an uninterrupted 
supply of their items. In recent years, the whispering problems 
of buyers being the non-availability of some durable products 
and their sub-components in time, Evehjthough industries have 
been using mathematical models, market survey methods to pre- 
dict the demand for their items in beforehand, there is no 

' 

single co st- effective model to give exact demand for a product. 

In an unstable world where industries are appearing and dis- 
appearing in a short period of time, the demand for a product 
of particular brand goes naturally to some other manufacturer, 
resulting in unpredictable demand fluctuations which are reported 
to be as high as 30 percent of the calculated demand in rare 
cases. 

In such a market environment fraught with uncertainties 
of demand and supply, substituting an item for another item with 
similar characteristics with inferior or superior qualities, 
becomes inevitable to maintain goodwill of the customer and to 
build the rapport. The above said environment typically i 

prevails in Indie® Market. 
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Eventhough. it is in this context that the problem of 
substitution is important, industries want to be cautious about 
the substitution policy. In an effort to improve goodwill of 
the customer and reputation, organisations do^aot want to take 
the risk of over substitution which has got its negative 
aspects also. By putting some policy restrictions and by 
maintaining an uninterrupted supply of the items in market, 
in the aforesaid situations, they avert problems. 

1.2 ORGANISATION OF THESIS: 

The present study deals with single item as well as 
multi -item inventory system under substitutable demand situations. 

In Chapter II, for complete substitution of two items 
for each other where items are of decaying in nature, a 
mathematical model is develops!, and solution methodology is 
given. A practical situation of a single item inventory system 
wlt&' emergency set-up cost, reworking cost is dealt and a 
simulation model for finding optimal inventory policy to meet 
a required service level is given in Chapter III. The Problem 

of N-item substitution with (Q,r f T) model is described in 

values of 

Chapter IV and solution methodology to find the/ decision variables 
the periodic review model with a simulation approach is explained 
in detail* 

For all the models dealt in this present work, numerical 
examples are given and results are analysed. 



CHAPTER II 


TWO ITEM SUBSTITUTION WITH BEGAY 

2.1 PROBLEM STATEMENT: 

Many organic compounds* medicinal items and. fashion goods 
deteriorate as time passes* This is reported to be due to the 
presence of some chemicals without which the compounds cannot 
be prepared. Sometimes the deterioration is accelerated by 
the amount stored and/or by other catalysts present in the items. 

We consider an inventory system &f a retailer of such 
items which are of decaying in nature. Among them* let us 
assume that few items have similar characteristics in terms 
of their utilities* purchase costs etc., and can be put in a 
group. We now concentrate on a grot?) of two items which can 
substitute each other. It is assumed that here in the situation 
of stockQput of one ites^its demands can be completely satisfied 
by the other. And there will be atleast one item in the stock 
at any time to avoid shortages. 

In this chapter, a simple inventory system with constant 
demand rate* infinite rate of production and variable rate of 
deterioration is considered where items are of substitutable 
in nature. A special fora of two parameter Weibull distribution 
is used to represent deterioration rate* which is reported to be 
realistic by Deb end Chaudhuri[2 ]. The demand of one item Is 
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completely fulfilled toy the other item and vice versa* In 
other words, at any time the system has the stock of only 
one item. 

A model is developed with the assumptions mentioned 
above and a numerical illustration is presented, 

2.2 LITERATURE REVIEW! 

Inventory models for items deteriorating with time have 
engaged attention of researchers in recent years. Various 
attempts have toeen madd to account for decay of items in 
inventory systems. An EOQ model for an inventory with a consr- 
tant rateJUertoratlon was developed by Ghare and Schroder^]. 
An order-level inventory model for Items with uniform rate of 
deterioration was introduced toy Shah and Jaiswal [6 3. In the 
recent publication, an EOQ model for Items with finite rate of 
production and variable rate of deterioration has been reported 
toy Deb and Chaudhuri [ 2 ]* Eventhough, many papers have toeen 
published, based on deterioration, literature is lacking in 
substitution of items which are deteriorating at a variable 
rate, 

2.3 ASSUMPTIONS* 

The following assumptions are made in designing the 

model. 

1) Demand rate of each item is deterministic and static. 

2) The replenishment is instantaneous . 
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3) At the time of procurement of item 1, stock levels 
of items 1 and 2 are zero. Similarly at the time of 
procurement of item 2, stock levels of items 1 and 2 
are zero. 

4) No shortage is allowed. As items are mutually inter- 
changeable complete substitution of one item by other 
leads to no shortage situation. 

5) Planning horizon is infinite. 

6) Unit variable cost of any items does not depend upon 
the replenishment quantities. In other words, no 
quantity discount is allowed. 

7) The substitution cost is directly proportional to the 
number of units to be substituted. 

8) Rate of deterioration follows a special form of Weibull 
distribution to make the problem mathematically tractable. 

2.4 NATURE OF DETERIORATIONS 

Two— parameter Weibull distribution is taken to represent 
deterioration rate, 

Weibull distribution function is given as, 

«{t) - t > 0 

where <^and p are positive constants. 


To make the problem mathematic ally soltfablSj a special 

form is considered by putting 0*2, we get, 

. ,2 

«(t) 


( 1 - -2M . 


) 
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By dropping higher powers of t (as (>t) n •* 0 asn**» anci 0) , 
we get, 

$(t) = at t >_ 0 

where 2 ^ = a. 

Here 0 < a < < 1 for practical purposes [2], 

Thus deterioration rate is linear with respect to time. 

2.5 NOMENCLATURE: 




Demand of item 1 
Demand of item 2. 

Deterioration rate of item 1 at time t t 0. 

Deterioration rate of item 2 at time t t > 0, 

Holding cost of iten. 1, in Rs./unit/unit time* 

% 

Holding cost of item 2 in Rs . /unit/ unit time* 
Unit cost of item 1. 

Unit cost of item 2* 

Flared ordering cost of item 1 in Rs*/order. 
Flared Irdering cost of item 2 in Rs. /order. 
Substitution cost of item 1 when it is 
substituted by item 2 in Rs./unit. 



Substitution cost of item 2 whoa it is 
substituted by item 1 in Rs./unit. 

Loss per unit item decay of item 1 in Rs./unit. 
Loss per unit item decay of item 2 in Rs./unit. 



A im I A I It 
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Ti me 


I-j-I nstantaneous level of inventory of item 1 
^“Instantaneous level of inventory of item 2 


Fig. 2*1 Two way substitution of two items 

(For items of decaying nature) 


[ 



8 


Decision Variables ? 

0^ The replenishment quantity of item 1 per cycle 

in units. 

Q 2 The replenishment quantity of item 2 per cycle 

in units. 



TC Total annual cost of both items consisting of 

inventory holding cost, procurement cost (fixed 
and variable) and substitution and loss of decay 
costs. 


2*6 FORMULATION* 


Fig. 2,1 shows the behaviour of on-hand inventory with 
time. At the time of procurement of item 1 stock level of both 
items are aero. Inventory level of item 1 is raised to 0^ and 
demand for both items are met with item 1* 

After stociQ 511 * item 1» item 2 is procured* This will 
meet the demand of items 1 and 2 for the period of duration T 2 » 

This cycle of length + T 2 is repeated subsequently. 

Referring to Fig. 2.1* the instantaneous level of 
\s 

inventory A derived as follows. Let O^Ct) be the inventory level 

at time t (0 < t < T.) * During the time interval (t* t+At) # 

' Ab 

the inventory level will deplete by an amount 9^(t),€i^(t)* due 

to decay and fey an amount (D-+D-) At due to demand. Thus, 

. Ab 

(t+At) m QjXt) ** ^(t).Q^(t) — (IX^+Dg) At 

Transposing and dividing by At, we obtain. 
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Q^t+At) - CL^(t) 
At# 


- ej(t) c^Ct) - 0\+i> 2 ) 


Taking limits as At ^ 0 yields, 
dCL (t) 

—|t:" * - V*) %. (t) ~ (%+V ° £ t < t- 

By rearranging we get. 


dCL (t) 

-If- + «iW <\ (t > 


(B 1+ D 2 ) 0 1 t < T ± (2,1) 


As the decay of item 2 starts only after its procurement 

the instantaneous level of inventory for item 2 is determined 

for 

following the same approach as/item 1 aid is given below* 


dQ 0 (t) 

-~lr~ + ° 2 (t) * 

Total time per cycle. 

By substituting ©^(t) « t 
and (2.2) we get, 

+ a^t Q^(t) « 

+ © 2 t ” 


dQ^t) 

~“ar~ 

dC^(t) 

IT 


- (D 1+ D 2 ) 0 < t < T 2 (2.2) 

T » t x + T 2 , 

and © 2 (t) » a 2 t in equns. (2.1) 



(2v3) 

- <vv 

(2.4) 


Solving the differential equation (2.3), we get. 


%(*) 


*w(a. ■+!>„) { 

4mm 



} [t + -gk t 5 ] 4- 



where i® constant of integration. 
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At t * 0, Q^(t) a» 0^, With this initial condition, 

» 0^, Substituting In eqn. (2.5), and substituting 

CX-t 

the expansion for exp (- -sp- t^), we get. 


0^(t) *s - (D^+D^) (t ~ -y t 3 ) + (1 - -y t 2 ) 

At t * T f G^(t) « 0. Substituting this in eqn. (2.5), we get, 

^ (D 1 +D 2 ) T 3 - ^ Q^ 3 - (Qj+B^) * 0 


After simplification, we get, 

Tj 3 + jiT^ 2 + + r » 0 


(2.6) 


where. 


p * -SOj/teCaj+Bg)) 
q * -3/a^ 
r ■* 30^/ (a^ ( D 2 +D 2 ^ ^ 


By substituting » xy ^ in eqn. (2# 6), we get. 





ax^ + b 


a* 0 


where. 


a 


§ (3q~p 2 

»»Q,^ + A(B» 4 -Drt) 

-3 ( -y 

%i 


b «* yr (2p 3 — 5pq + 27r) 

g^Cl^^) 2 - q^ 5 

Ao^ (E^^D 2 ) 3 


(2.7) 
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0 is 


Solution for + ax^ + b 

A_ +B n A» “B_ r — A-. 4B-| 

Xj, » VV * + “^2^ W~3) , — 


^(xpi) 

( 2 . 8 ) 


•where i 


*1 


B, 


/ b A 
V W + 


/ b 

C 2“J 


b 2 

T 27 ' 


b 2 ^ 

; 


After simplification, we get. 


^ » [x^ 5 ^ + jK 3 a L 6 ^ 4 Q 3 L 4 ^ 5 % 2 + K 6 1 

B 1 “ ^' K 1 Q 1 5+K 2 a i * s 3 


1/3 


1/3 


(2.9) 

( 2 . 10 ) 


where* 

K-j^ *■ -1/8 (D^^^) 

Kg * 5/4 O^+Bg) 

Kj m 0.0 

K^ « *"3/(®(^"*‘^23 a i3 

- -3/ fr&ij 2 C 0 !^) 2 ) 

Kg * —l/ 8^ 

For making the analysis easy, rewriting Hie roots of 
the equation (2*7) in different form, we get. 


** y^*®^ — z^(Aj— B^) 
where, (y^, Zj) «an take 

t ** 

(l. 0,0.0) ©r (—0*5,0.866) or (~0.5 ,-0,866) 


( 2 , 11 ) 





Let A 11 » K^ 3 + (2.12) 

B 1X * C-lO^Oj 4 + K^Q^ 2 + Kg)) 1 / 2 (2.13) 

As constants to Kg are negative and values within 
square root terms of A^ and 

getting real roots after separating real and imaginary parts. 


B^ are negative, «e will be 


Rewriting values of and in eqn. (2.9) and (2.10) 

we get, 

\ m (A^ + B^) 1 / 6 (Cos (|) + i Sin (|)) 

B 1 “ (A 11 + B ll^ 6 ^ Cos “ i Sin (|)) 
where, C « tan"* 1 (B^-j/A^) 


Finding values for A^+B^ and A^-B^, we get, 

A 1 +B 1 - 2(A 2 1 +B 2 1 ) 1 / 6 Cos (C/3) (2.14) 

A 1 -B 1 » 2(A 2 1 +b| 1 ) 1 / 6 Sin (C/3) (2.15) 

E <gu (2.11) will become, 

** 2y^ ( A^i +B 2 .j ) 1 /^ Cos (C/3) 

- 2^4(4^^ ) 1/6 Sin (C/3) (2*16) 

By putting we can find 3 values of T^. 

For item 2, solution methodology for solving the eqn. (2.2) 


is same as for eqn. (2.1). 



finally we get. 



lafeere , 

B 21 

C 1 

*2 

®3 

R 4 

% 

% 


\®2 + R 2 Q 2 
(-(R^Qg^ + R^Q 2 2 + R g)) 


1/2 


-1 


B, 


tan ^ v^ 1 
A 21 

- i / sCi \+ d 2 ) 

3 / 4 ( nj + D 2 ) 


«= 0*0 

■ a 2 
*» *-3/ (16(0^ +0 2 ) a 2 ^ 

** * 1 / ttg 


and (yg*^) P®* 1 ' fcake value# (l.0 t 0.0) or (-0 

(**0 • 5 ,*• 0 * 866 ) *> 

A CL 

By substituting “ x 2 ~ 


ve get 3 values for T 2 * 



Formulation of Total Cost: 


(a) Procurement Costs 


For item 1 » 


For item 2 


A 2 + C 2 Q 2 


(b) Inventory Holding Costs 

T- 

For item 1 » 


I hi^Ct) 


dt 


0 


/^[Q^CL— -ip t 2 ) - (Z^+E^Xt- -ip t 3 )] dt 


0 *^ 2 4» ^ Ct*J /. . 

[C^Ct - t 3 ) - (D 1 +D 2 ) ~ t 4 )] 


“A 


T a 2 


»i c°i( T i - - < WH - ^ir )] 


For item 2* inventory holding cost 


aJ£* T 2 

h 2 lQ 2 (T 2 - 


(c) Substitution Cost* 


For item 2 
For item 1 


V„ D„ T- 
2 2 1 


V 1 “l T 2 
(d) Loss due to decay of item 1 


Sit 


f *F l0l [<^(1- ^5-) - dt 

F l“l I l CC b. (1 " "V 1- ) ~ 
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For item 2, loss due to decay, 

~ F 2« 2 T 2 £ q 2 &~ ( D i +D 2 ) ( “T “ 

Thus Total annual cost 

a I 5 

TC " + kjA (T i “ “V 1- } 

a 3 3j2 a y 4 

* h 2^2 < S 2 * -¥-> ~ WV ( V “ -fer^ 

- V W <%- - ^-) - w 2 + Wi 

a T 2 ^ j 2 

* F M^ a - -V- > + F 2“2 Q 2 I 2 & - “V-J 

T. a-,*, 5 

*■* ^2 a l^I ^l + ^2^ ^~2* ** ^ 

T g j 5 

- F 2 a 2 T 2 (nj+Dg) (-£- (2.18) 


Setting the partial derivatives of TC with respect to 
and Qg equal to zero yields. 



and. 


(▼•& 5 T ** u m 0 


IS^ ■ (t I^* u ^ )/vZ ■ 0 


where, 


W at 


Tj+T 2 


u *• Aj^eAg + *i* CgQg + (f^ 

3 


g A \ 


(2,19) 


( 2 , 20 ) 


+ ^°2 (t 2 - ~ s r~ ) * W*V Hr* -^g-) 

*P 2 a T ^ 

- j>2 cn^) (-4- - -tJ-) * w 2 
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Bu 



a 2 T l 


+ ^2 B 2 T 1 * F l a lVi ^ IT~^ 
+ ( 1- ^ ) 


•k* c i*x 


T- 3 


F l a l^l A +B 2^ ^ 2 12 ^ 

T 0 . a^Z 3 

F 2 a 2 2 '2 ^ D l +0 2^ A T?"'^ 


lA A 


cl+; 


a 


A 


3cti^ T, 


tr~) + *iA (1 - -p-> 


4a n T. 3 


- h, <vV (t, - -ir-) 

T a, TL 3 

+ ? 2 0 2 + F i a i A * Ar-* 


+ F x a A (1 


*lV *- 


) ~sr ■ 

F l a l A +1 V A - — 


a 0 T 0 ^ 3a «5^ 


c 2 + ^ (t 2 ■ *^r~ ) + Vfe } ~i£ 

4a l 1 3 5 T 

- »2 <w <*2 - - fr-j -<4 + 7 A -*§ 

* F 2“2*2 & ~ ^ * F 2“2®2 <* * -S 


-m 5 3 5? 

F 2®2 A* 1 ^ A * HEr~> 


( 2 * 21 ) 


(2*22) 
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d(T 1+ T 2 ) 

6% 


8? 


(2.23) 


c)v 

set 


a(T 1+ T 2 ) 


(2.24) 


where, 


i . __ <v®i> _ ( W ^ 

3~?tr 2 T + y i Sf H *5[ U * 25J 


^ - 2 C I <4l + B^)-^ 


(2A n (3KiO l 2 +K 2 ) - (4K 4 0 L 3 +2K 5 O l )) 

Cos (^) + ( (~Sin (^) ) ^ 

“ ®UL "^$£^*11^ 


( 2 . 26 ) 


•~^ g r ~ “ 2[ I (A ia^i )_5/6 


(2Au<3Ki Q 1 2 -ric 2 ) - (4s J( a l 3 4aE 5 a L )) 
Sin (^) + | (^Ii+bIi) 1 / 6 Co® (^) 
(Cl*(^) C&n - ^ 


(2.27) 
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where, 


d. 


dB ll 

C 


5Kl°l 2 + K 2 

(2.28) 

-(AK^ 3 + 2K 5 -Q 2 )/(2B 21 ) 

(2.29) 

n B 11 
tali’ (f”) 

*11 

(2.30) 


Similarly 

dT 


1 (A 21 +B 21 ) 

* y 2 — s- — * z 2 


^ A 21 ~ B 21 ^ 

(2,31) 


where. 


- 2 1 1 <4. ♦ B ii r5/6 

(£A 2 i (3it L Q 2 2 + Rg) * (AR^Qg 3 + 2HjC^)) 
Cos ( ^ ) ♦ (A^ ♦ B ^) 1 / 6 (-Sin 


(( 1 + (j^) 2 )*"^) ((Ag^ — | 


B 



21 — SCT>/A 21 )] 


(2.32) 


where. 




dB, 

dL 

c i 


31 ^ Qg 2 ♦ Rg 

**(4R^. Qg^ + ^Q;>)/ (2By| ) 

taET 1 (^ t ) 

*21 


(2.33) 


( 2 . 5*0 


(2.35) 
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Using eqns. (2.21) - (2.35), and substituting the terms 
in eqns. (2.19) and (2.20) we get. 


bTC 


m 0 


dTC 


*2<W 

0 


(2,36) 


(2.37) 


Solving the non-linear equations (2.36) and (2.37) by 
simultaneous iterative method for various values of and Qg 
we get various combinations of Q^, G^, and T 2# which will 
give local optimum solutions. 2he values of €^,0^, 2^ and T 2 
which give global optimum are found. 


One problem is solved and the results are tabulated. 

Exhaustive enumeration method is followed, because: 

(i) Two local minimum points might be giving same TC values 

at 

but one might give/low values of &. # Q p ,T 1 and T ? 
comparing to the other. In this, manufacturers will give 
preference to the point which comes at low values of 
decision variables. 

(ii) Many a time optimal point will not be at * 0 and/or 

~~~ o 0. So only scanning the whole function will be the 
best approach. 

(Hi) Problem in using non-linear programming techniques is 

that finding initial points to start iteration Is found 
to be difficult. 
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Table 2. It Total cost vs. 0^ and Qg. 


5 

°2 

Total 

Cost 


T_ 

2 

130.00 

10050.00 

5728.48 

26.67 

60.07 

150.00 

10050.00 

2891.82 

26.563 

60.07 

170.00 

10050.00 

3938.75 

26.46 

60,07 

190.00 

10050.00 

4989.10 

26.35 

60.07 

210.00 

10050.00 

6042.82 

26.24 

60.07 

230.00 

10050.00 

7099.83 

26.13 

60.07 


Total ©oat la minimum at » 150, * 10050, 

T x * 26.7, T 2 * 60.07. 



CHAPTER III 


SIMULATION MODEL FOR ONE WAY SUBSTITUTION 

3.1 NATURE OF THE PROBLEM* 

Let us consider a manufacturing situation where end- 
item is In a fluctuating demand. Mathematical models show that 
there could be an error as high as 30 percent in the forecast 
of the demand for the item. We assume that the product can 
be substituted by some other available item with similar 
characte ri st ic s. 

Assume that instead of buying it at the moment the demand 

exceeds inventory level* the industry plans to keep an inventory 

of the product which is going to substitute the original. To 

tackle the problem of non-availability of space and to achieve 

some economy by making use of the difference in cost of First 

s 

level items (end "items to eustomez^ and second level items (items 
which are one stage prior to the end -items ), j^ssiase that the plan 
is to buy the end-item substitute and/or second level items. So 
the company has to rework on second level items before assembling 
it to make the end-item substitute, if second level items are 
also purchased. 

The selection of buying end-item substitute and/or 
second level items depends .on the cost of end— items# end— item 
substitute, second level items, holding cost of first and second 
level items, emergency set-up costs, and cost of reworking and 
assembling end-item substitute from second level items. 
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Our objective is to meet a particular service level by 
substitution at an economical inventory carrying cost. 

3.2 LITERATURE REVIEW? 

Many attempts have been made by researchers to solve the 
problem of unsatisfied demand by substituting it with an item 
of similar characteristics. In a paper by Carlson and Yano [ l], 
the formulation of an algorithm for fixing buffer levels taking 
into account the uncertainty of demand in MRP systems is being 
quoted. Eventhough there are many mathematical models for 
substitution of demand developed by Nigs© C 5l» easy-to-use, 
simple and cost-effective methods are lacking in this area. 

In practical situations, it is very difficult to represent 
a system with many constraints putforth by market conditions and 
organisation policies mathematically, because of the analytical 
complexity involved in modelling. With this in view, in this 
chapter a simulation model is developed to find the optimal 
inventory policy of an inventory system taking into account the 
factors like emergency set-up cost, reworking cost and inventory 
holding cost of end-item substitute and second level items at a 
specified service level* Numerical examples are given and 
results analysed* 

3.3 SYSTEM DESCRIPTION? 

Mm there are three different situations assumed, we 
consider here three different ways of working of the system. 

They are* 

(l) Situation where reworking to make end-item and emergency 
ordering of end— item s proves to be economical, In this situation 



the required, quantity of end-items and second level items are 
procured. The customer's demand for item considered is fulfilled 
with the same, and side by side, the end- item substitute will be 
in the making from second level items# If the desand for end- 
items exceeds, it will be fulfilled by the substitute which 
Is made available by reworking and assembling. But if the demand 
goes beyond the procured quantity, we go for emergency ordering 
of end- item substitute. If enough time is left, emergency 
set-up is carried out till fill-rate crosses the required 
service level. 

In case of actual demand lesser than procured quantity, 
the excess is carried over to next period and that quantity 
of excess items will be reduced in the quantity to be ordered 
for third period* 

This procedure is followed subsequently. 

(2) Situation where stocking of end-item substitute and its 
emergency ordering proves economical. 

In this situation only end-items and end-item substitutes 
are procured in inventory. Initially the customer demand is 
fulfilled by the end-item and then by end-item substitute. If 
the demand goes excessively, emergency set-ups are made. 

In case of demand lesser than inventory, excessive end- 
items and end-item substitutes are carried over to next period 
and the ordering quantity for third period Is accordingly 
adjusted. This follows for subsequent cycles too. 



(3) Situation Viter© stocking of end— item and end-item substi- 
tute proves to be economical i 

Demand is fulfilled by end-item , end— item substitute 
nventory carrying and/or emergency set-ups are as in 
cases 1 and 2. 

3.4 ASSUMPTIONS! 

Following assumptions are made to characterise the model. 

1. One end- item is made of two second level items each. 

2* Demand forecast for endr-itera in all periods of the 

planning horizon is available. 

3* Demand for end-item are in critical situation. 

4. Individual demands for second level items are not 
considered. 

5. The forecast of demand is equal to the mean demand and 
forecast errors are assumed to be distributed normally 
with mean as zero. 

6# The standard deviation of the forecast errors for an item 
in a period is assumed to be a constant fraction of the 
mean demand of that, item in that period. 

7. The holding costs, for a period are charged on the inventory 
carried oyer from that period to next. 

8. Mo other uncertainty exists in the system* lead times and 
quantity supplied are assumed to be detemLnlsrtic* 

9. Emergency set-ups are allowed only to procure end-level 
items. 

10. Lead t Sxm for emergency set-ups are finite and lesser than 
the lead time in case of regular set-ups. 

1*1 r -et-un cost is higher than set-up costs per order 



which are encounteved in regular procurements. 

3.5 SERVICE LEVEL MID STANDARD DEVIATION; 

There are tiro principal approaches to the selection 
of safety stock, namely specification of a desired service 
level and explicit costing of shortages. Two popular measures 
of the former aret* 

(!) Specified fraction of demand to be satisfied routinely 
from stock which is also called the fill-rate. 

(ii) Specified probability of no-stockout per replenishment 
cycle. 

Vie will be using the former one because of its practicality 
and easy-to-use nature. 

As forecast errors are very difficult to estimate even 
under stationary or slowly ©hanging situations* we assume that 
the standard deviation of forecast errors is a constant fraction 
of the demand forecast which is the mean of the distribution. 

This has been reported to serve well in many situations by 
Silver £ 7]* To make the problem tractable a relatively simple 
method to account for the changing nature of demand and hence 
the assumption* 

3.6 NOTATIONS! 

C - Fraction of demand giving the standard deviation. 

T - H as her of periods in the planning horizon. 

K|l) - Safety factor for item i* to achieve desifed service 
level there i « 1 for end-item substitute, 

i * 2 or 3 for second level items. 
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Mean of the demand for item i in period J. 

Standard deviation of demand for item i in period J. 
Safety stock required for item i in period 
Order upto level for item i in period 3. 

Major set-up Cost for end- item. 

Minor set-up cost for end-item. 

Minor set-up cost for end-item substitute (i**l), 
second level items (i «* 2,3). 

Inventory holding cost for one unit of end-item. 
Inventory holding cost for one unit of end- item 
substitute (i»l) , second level items (i-2,3). 
Quantity of item i to be ordered in period, 

Quantity of end-item i on hand in period j. 

Cost of reworking involved in making end- item from 


second level items in Rs./unit of end-item substitute. 
Emergency set— up cost in ordering end-item substitute 
in case of emergency in Rs./N items* 


An emergency order for H items or below, emergency 
set-up cost will be one EMRGEN. For an emergency order 
for X items where N < X < 2N, the emergency set-up cost 
be to g two EMRGEN and so on* 


mmm 
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Fig. 3*1 (a)* BLOCK A OF FIG* 3.1 
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3.7 HEURISTIC ALGORITHM: 

Joint lot sizing algorithm/ a simple version of lot- 
sizing algorithm (Vagner and Whit in [8 ])) which was developed 
by Gopinathan [4 ] is used along with Dynamic Rolling Schedule. 
The functioning of Dynamic Rolling Schedule is given in system 
description , whereas the heuristic algorithm is given in 
Appendix A* 

3.8 SIMULATION MODEL: 

A Simulation Program is written in FORTRAN with subroutirn 
HEURIS and GENER, to find the optimal inventory policy of how 
much to procure and which type of items to be procured to meet 
the required service level. Fig. (j .1) shows the flow chart of 
simulation model. 

HEURIS* 

This subroutine finds quantity to be procured in terms of 
end-item , end-item substitute and/or second level items by 
applying the heuristics mentioned earlier in this chapter. This 
Is done before applying Dynamic Rolling Schedule. 

GENER* 

This subroutine generates Random Numbers (which are 
distributed normally) to find actual demand for the given mean 
demand and standard deviations. A system supplied Random Humber 
Generator G$5fj3DF from NAG is used by this subroutine. 
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Objective Functiona l 

Objective is to find the optimal ordering policy* 
given the reworking costs* emergency set-up cost, major and 
minor set-up costs for a required service level. 

3.9 SIMULATION EXEEUMENT: 


"Problems are solved using the simulation model 
described above with different cost parameters for an 
inventory system of 

Planning Horizon 1 year 

Number of periods 12 

Emergency set-up time l/ 5 th of a month 

Eili-rate 90 percent of demand. 


The input parameters for the ' problems are given 
in Table 3*1, and 3.2 and the output results are shown in 
Table 3.3 and 3.4. 

From the output results the following observations are 

made. 

1 ) As the reworking cost is low* it is obvious that better 
to maintain second level item inventory instead of going for 
end— item or end-item substitute inventory* and emergency set-ups. 
The graph (Fig. 3.3) shows that the effect of K(l) and K( 2 ) on 
total cost* at 90 percent fill rate. 

2 ) For the cost structure given in problem 2 * at zero level 
safety stock for all items, it is found that emergency set-up is 
economical* to maintain a fill-rate of 90 percent. Fig. 3,4 shows 
at zero level safety stock of second level items, and with an 
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increase in safety stock of end- item inventory, there is a 
decrease in total annual cost. The trade-off "between emergency 
set-up cost and inventory holding cost is at K(l) » 5 percent and 
after that the total annual cost started rising with increase 
in safety stock level of end-item substitute. 
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Table 3.1 1 Input values for the Problem. 


Cost terms 

Rs. 

flavor set-up cost 

50.00 

Minor Set-up Cost for 

End-item 

25.00 

End- item Substitute 

25.00 

Second level item 1 

15.00 

Second level item 2 

10,00 

Holding Cost for 

End- Item 

1.10 

End- item Substitute 

1.10 

Second level item 1 

0.53 

Second level item 2 

0.40 


Table 3.2* Input values for the problem. 


Problem Mo. 

Reworking Cost 

Emergency Cost 

1 

0.6 

1.4 

2 

1.4 

1,4 
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Table 3.3* Total cost vs. K(l) and K(2) for Problem 1. 
(In Rs.) . 


K(2) 

K(l) 


Total Cost (Rs. ) 

0 

5 

10 

0 

1882*28 

2748.63 

3388.55 

5 

1385.18 

2791.53 

3431.45 

10 

1425.13 

2829.48 

3469.40 

15 

1527.68 

2874.03 

3513.45 

20 

1626.73 

2913.08 

3553.00 

25 

1729.63 

2955.98 

- 
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Table 3*4? Total cost vs. K(l) and K(2) Tor 
Problem 2 (in Rs.). 


K(2) 

K(l) 

0 

5 

10 

Total cost (Rs.) 


0 


1665.75 

1702.23 

1733.95 

5 


1577,10 

1613.67 

1645.30 

10 


1607.03 

1643.50 

1675,22 

15 


1637.42 

1673.95 

1705.63 

20 


1665.93 

1702.33 

1734.13 

25 


1697.27 

1733.67 

- 

30 


1729.20 

1763.50 
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Fig. 3-5 Effect of K(1) and K(2) on total cost 



CHAPTER IT 


SIMULATION FOR H-ITEM SUBSTITUTION 

4,1 PROBLEM DEFINITIONS 

We consider an inventory system of N-items, We further 
assume that among these items some items are of interchangeable 
in nature* It means, in case of stock-out of one item, its 
demand can be fulfilled with some other item(s) of similar 
characteristics. If no such item is available, its demand is 
backlogged and fulfilled in the next period* 

In this chapter, a simulation model is developed for 
widely used periodic review model, generally called (Q,r,f) 
model. As this model precludes an analytical approach, with 
substitution, simulation technique seems to be a suitable 
method in evaluating the decision variables of the model. 

The objective is to find the values of Q,r* R and T of the 
(Q,r f T) model at minimum total cost with N-item substitution for 
the given set of cost parameters* 

4*2 (CLr.Ti HPBKH 

The decision variables of the model ares 
Q — Quantity to be ordered whenever an order is placed. 

r — Reorder level — the quantity which determines the need 

for ordering or not ordering. 
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T — Cycle time - a stipulated period after which the 
status of inventory is reviewed. 

After a time period T, inventory is reviewed and if the 
inventory position is below reorder level, an order for quantity 
Q is placed. 

The (Q,r,T) model for N-item with two-way substitution 
can be explained as follows, \*henever an item I D is demanded, 
it is satisfied if it is in stock, otherwise it is satisfied 
by any item Ig, the most economical of the items that can be 
used for substituting where S #» l,2,.,.,n and B /£ S»from set 
of n items. If the item I D *s demand cannot be satisfied fully 
by itself or by any other item, then the remaining quantity 
is backlogged. : 

4.3 ASSUMPTIONS* 

Following assumptions are made in characterising the models | 

1) Only certain items can be substituted for other item of 
similar characteristics. 

2) The unsatisfied demand for an item is fulfilled by subs- 
titution if possible. The unfulfilled demand is backlogged. 

3) Leadtime for procurement is static and known# j 

4) Parameters of the distribution functions are know. 

5) Demands are stochastic and follow geometric distributions. ! 

6) Inter-arrival time for item - demand is exponential . . j 

7) Planning horizon is known j 

f 

Other aspects considered are discussed in the following j 

# l 

sections. ' 
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Can-Order Level ? 

Mien the inventory is reviewed, if the inventory position 
is near to Re-order level hut not below it, according to (Q,r»T) 
conventional models, an order is not placed for it. It Is 
obvious that in a N-item inventory system, if order is not 
placed for it, it may increase the substitution cost as well 
as backlog cost. Hence can-order level is also considered* 

At the time of review, if the inventory position is below can- 
order level and/or above reorder level, an order is placed. 

Priority Order s 

Three main events in the model ares 

1) Customer Arrival 

2) Arrival of Procurement 

3) Review 

Mien two or more events occur simultaneously, the events 
are selected according to the following priority rules, 

1) Arrival of Procurement 

2) Review 

3) Customer Arrival, 

Sollt-up of Coat s 

For all items, the split~up of cost are given as, 

1) Ordering cost. 

2) Holding cost. 

3) Backlog cost. 

4) Substitution cost if applicable. 
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4.4 NOTATIONS* 

The notations used in developing the simulation model 


are* 



FQTY(I) 

mm 

Fixed ordering quantity for item I. 

TIM 

- 

Time at any point. 

T 

- 

Cycle time. 

REOD(I) 

- 

Reorder level for item I. 

CNOR(I) 

* 

Can-order-level for item I. 

HOLD (I) 

mm 

Holding cost of item I. 

UNTC(I) 

- 

Unit cost of item I. 

BCLG(I) 

- 

Backlog cost of item I* 

SUB(I,J) 

— 

Substitution cost if item I is substituted by 

some other item. 

LEDTIM(I) 

- 

Lead time for item I. 

TC(I) 

mm 

Total cost involved in maintaining inventory of item I. 

INV(U) 

- 

Inventory level of item I» in period J. 

DMD(I t J) 

- 

Demand for item I in period J . 

TLBfT 

* 

nT-TIM. Time left in that particular cycle h. 

LREOD 

mm 

Status of Reorder level. 

If LREOD « 0, Order is yet to be placed for next period 


LREOD > 0, Order is already placed, 

PROMs - Parameter P of demand distribution. 

AKVR - Parameter of inter-arrival distribution* 

4.5 METHODOLOGY* 

The flow chart of the simulation program is shown in 
Fig. 4.1, 4.1(a) and (b). 

























FIG. 4.1(b): BLOCK D OF FIG. 4.1(a) 




To start with, the EOQ for each item is found assuming 
that all are ordered separately and the various time periods are 
calculated. This EOQ is called fixed ordering quantity, and 

g o -Q 

denoted by Q . 10 percent of the Q are fixed as r (reorder 

s s 

level of item S) and 30 percent of Q are fixed as R (can-order 
level of item S). 

Among the set of cycle times calculated for various items, 
minimum cycle time is chosen, and simulation is done for item I 
for the calculated Qjr,®- values as mentioned earlier. 

Simulation is actually generating time values and corres- 
pondingly demand values for Item I. Time (TIM) is check®! to 
determine the occurrence of the events in the model. 

Reordering of item 1 is done 

if the inventory level is below can-order level and/or 

if the time left in the cycle is only 3ust enough to 

procure item 1 for next cycle in time reordered now. 

If the inventory level comes below zero in a cycle due to 
fluctuations in demand, the item will be substituted with some 
other item(s) • The selection of item is based on availability 
end economy in the cost of substitution. The unfulfilled demand 
is backlogged said backlogging cost is also taken into account. 
This is repeated in subsequent cycles for item 1* Total cost of 
maintaining inventory is calculated. 

This procedure is repeated for all items with same cycle 
time. Total cost of inventorying all items is noted alongwith 
cycle time. 
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The procedure mentioned above is repeated taking different 
possible cycle times# and the total cost values are calculated 
in each case* 

£>^ock t- ( S 'Tor 'fbUowin^, 

%e cycle time which gives minimum total cost for 

inventorying all items is taken. Then cycle time is perturbed 

on both the sides by 10 percent and the simulation is repeated 

to find the time period which is giving optimum total cost. 

Now with the above initial parameters, the simulation 
is done for a planning horizon by reducing r^ by 25 percent and 
then cost of inventory policy is noted each time by increasing 
r^ by 1 percent, until r^ is 1.25 times the initial r-^* The 

corresponding to the minimum cost is fixed as the parameter r^. 
Then the simulation is done in the same way to find r 2 , T^ t **,r n 
values. 

The same methodology is adopted to find and 

® 1 ****» 

4.6 NUMERICAL EXAMPLE* 

A simulation program is written in FOKTRAH to validate 
the model developed. A. numerical example is solved and results 
are shown. 



NUMBER OF ITEMS — 3 

planning horizon — 5.00 

1 ’ 1 ’ 

ITEM NO . 

I FIXED COST f- 
35.00 

hold, cost 1 
4.5C 

BACKL03 COS T 
2.0o 

2 

30.03 

3.20 

2.80 

3 

40.00 

4.30 

2.P0 

SUBSTITUTION MATRIX IS 

AS FOLLOWS 


ITEM NO |, 

‘ , • 2 

3 

, 

H 

99. 99! , 0.901 

0.85! 


21 

, 99.991 

0.75! 


31 

O.sol 0.80 1 

9°. 991 


DISTPIBUTION details 



ITEM No f 

L£4DTT«: .1 DEMAND 1 INT-ARV 

-TIME | 

11 

0 .1 0 f 0.921 

0.50 1 


21 

O.ior 0.951 

0.701 


31 

O.lOf 0.991 

0.651 



10.00 

9.00 

12.03 


*** END 3 P INPUT DETAILS *** 


OPTIMAL PARAMETERS OF THE POLICY WITHOUT SUBSTITUTION 


OPTIMA!. REVIEW PERIOD -- 0.54 YEAR 

PLANNING HORIZON IS 5.00 TEARS 
"UMBER OE ITEMS -« 3 

following are optimal ouantities 

ORDERING 3 'JANrirlES 


ITEM NO I OUANTTTlf 1 




i i 

2f>.05 f 

2 1 

25.13 r 

3 f 

42.83 ! 

CkH 0P3 

EP LEVELS 

ITE.« mo ) 

DU A N TIT! J 

1 t 

5.45 f 

2 1 

1.24 1 

3 l 

12.47 1 

REORDER LEVELS 

ITEM MO 1 

Du^sriTr r 


1 1 

2.9a 1 

2 1 

2 . °2 « 

3 I 

4.68 f 

FOLLOWING ARE THF COST ESTIMATES OF TrfF. POLlCr 

TOTAL COST 

IS 1702.14 


cosT-s»ijr-yp is as follows 


ORDERING COST § 50.00 
HOLDING COS’- J§|-f I 
BACKLOG COST 555.48 



CHAPTER ? 


CONCLUSION AND SCOPE FOR FUTURE M}RK 

5.1 CONCLUSIONS* 

In the present study, various deterministic as well as 
stochastic models have been developed for two-item/ multi-item 
inventory system to take account of Interaction of demand due 
to shortages and/or introduction of new items. The models 
take account of many realistic situations in inventory manage- 
ment. Some numerical examples have been taken to examine the 
possibility of substitution. 

For deterministic demand case, the model with independent 
procurement of two items has been developed. The complete 
substitution of one by the other in case of stock-out of the 
former is discussed. It is found that the policy is to go for 
a very lengthy planning horizon* A simulation model is developed 
to find the optimal ordering policy for a single item case. 

In this model, the possibility of substituting one item by 
other item in one way is discussed and solution methodology is 
given to determine the policy of whether to go for end-item 
substitute or to bfy second level items and reworking to make 
end-item, is given. It is found that the optimum ordering policy 
is determined not only by emergency set-up cost end reworking 
cost but also by service-level. 



'■ : 

A simulation model for B-item with two-way substitution 
is studied and the steps involved in finding the values of 
decision variables which determine the optimum inventory policy 
are given. It is found that, to find global optimum , several 
simulation runs have to be made. Substitution proved to be 
economical sometimes, and costlier in some other cases depending 
upon the cost of substitution. 

5.2 SCOPE FOR FUTURE WORKS 

In order to account for realistic situation, models 
developed in this thesis need to be extended to account for various 
practical situations. 

Model developed in Chapter II, where two-item substitution 
with decay is considered , can be extended for items with different 
decay patterns* To be more realistic, quantity discounts can also 
be included. 

The model for single item, one-way substitution can be 
extended to N-item substitution to save storage spa©2a».d cost of 
maintaining by keeping second level Items for the end-items 
considered. In practical cases, there will be independent demand 
for second level items. So models accounting for that too will 
be of more use in practice. 
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APPENDIX A 


SIMPLE VERSION OF JOINT LOT SIZING ALGORITHM 

Initially lot-by-lot approach is assumed* i.e* the items 
needed for a particular period Is purchased at the beginning of 
the period. 

Now, starting from the last period, the trade-off between 
holding cost and procurement cost is made assuming that if the 
quantity needed for lastperiod is ordered in the penultimate 
period. If there Is a saving In that, the Items needed in last 
period are purchased in penultimate period. If there is no cost 
saving in procuring it in previous period, lot— by-lot approach is 
assumed for last period. 

This is applied to all periods starting from end period, 
and economic order quantities are determined* 
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